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ABSTRACT 

We report the results from our analysis of Suzaku XIS (0.5-10 keV) and 
HXD/PIN (15-40 keV) observations of five well-known local ULIRGs: IRAS 
F05189-2524, IRAS F08572+3915, Mrk 273, PKS 1345+12, and Arp 220. The 
XIS observations of F05189-2524 and Mrk 273 reveal strong iron lines consistent 
with Fe Ka and changes in spectral shapes with respect to previous Chandra and 
XMM-Newton observations. Mrk 273 is also detected by the HXD/PIN at ~ 1.8- 
o. For F05189-2524, modeling of the data from the different epochs suggests 
that the change in spectral shape is likely due to the central source switching 
off, leaving behind a residual reflection spectrum, or an increase in the absorbing 
column. An increase in the covering fraction of the absorber can describe the 
spectral variations seen in Mrk 273, although a reduction in the intrinsic AGN 
luminosity cannot be formally ruled out. The Suzaku spectra of Mrk 273 are well 
fit by a ~94% covering fraction model with a column density of ~10 24 cm" 2 . 
The absorption-corrected log[L 2 _io hev/Lm] ratio is consistent with those found 
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in PG Quasars. The change in the spectral shape on a time scale of a few years 
implies that the absorbing matter must be near the AGN (~1 pc). The 0.5- 
10 keV spectrum of PKS 1345+12 and Arp 220 seem unchanged from previous 
observations and their hard X-ray emission is not convincingly detected by the 
HXD/PIN. The large column density derived from CO observations and the large 
equivalent width of an ionized Fe line in Arp 220 can be reconciled by an ionized 
reflection model. F08572+3915 is undetected in both the XIS and HXD/PIN, 
but the analysis of unpublished Chandra data provides a new measurement at 
low energies. 

Subject headings: galaxies: active — galaxies: individual: IRAS F05189-2524, 
IRAS F08572+3915, Mrk 273, PKS 1345+12, Arp 220 — galaxies: starburst - 
X-rays: galaxies 



1. Introduction 

Ultraluminous Infrared Galaxies (ULIRGs) are galaxies with Lir, = Ls-ioooAtm > 1O 12 L , 
equivalent to the minimum bolometric luminosity of QSOs. The discovery of a large popu- 
lation of local ULIRGs is one of the most important legacies of the IRAS satellite. At lumi- 
nosities above 10 12 L Q , the space density of ULIRGs in the local universe is greater than that 
of optically selected quasars with similar bolometric luminosity by a factor of ~1.5 (Sanders 
& Mirabel 1996). Thus ULIRGs represent the most common type of galaxy at these high 
luminosities. Systematic ground-based and space-based optical and near-infrared imaging 
studies have shown that local ULIRGs are almost always undergoing mergers (e.g., Sanders 
et al. 1988; Veilleux et al. 2002, 2006). Sanders et al. (1988) suggested that these objects 
represent a dust-enshrouded phase that eventually evolves into optically-selected quasars. 
If this is true, ULIRGs take on a fundamental importance for the origin and evolution of 
quasars. 

The primary energy source — active galactic nucleus (AGN) versus starburst activity 
- of ULIRGs is still under debate. Optical and infrared observations show t hat most local 



ULIRG s are dominated by starbursts, but about 30% show evidence of AGNs ( IVeilleux et al. 



1999al lbh. Furthermore, the "warm" infrared colors (/ 2 5 ^m/ fm (im > 0.2) and quasar-like 



spectra of the more luminous objects imply that active nuclei are significant in these objects 
(e.g., Veilleux et al. 1995, 1999a,b, 2008; Genzel et al. 1998; Tran et al. 2001; Armus 
et al. 2007). These results would be definitive if it were not for the possibility of severe 
obscuration. 
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Galaxy mergers cause massive in-falls of gaseous material towards the center of ULIRGs 
(e.g., Rupke, Veilleux, & Baker 2008), and column densities > 10 24 cm -2 have been deduced 
from CO measurements (e.g., Downes & Solomon 1998; Evans et al. 2002). Observations 
at UV, optical, near-infrared, and even mid-infrared wavelengths may therefore not always 
probe the cores of these objects. High resolution (<^C 1"), high frequency radio observations 
can penetrate high columns and are an excellent probe of whether an AGN is present (e.g., 
Lonsdale, Smith, & Lonsdale 1993). However, the bolometric luminosity in the radio band 
is insignificant, so radio data cannot test whether accretion onto a supermassive black hole 
(SMBH) is the dominant energetic process. We are left with hard X-rays. 

Recent surveys with XMM-Newton and Chandra (e.g. Ptak et al. 2003; Franceschini 
et al. 2003; Teng et al. 2005) have found that ~40% of observed ULIRGs show signatures 
of AGNs. The observed 2-10 keV luminosity of the surveyed sample is ~10 40 -10 43 ergs s , 
with a majority of the sources having luminosities below 10 42 ergs s _1 . While the ratio 
log[L 2 _io kev/^m] is small in nearby ULIRGs (from —4 to —1; e.g., Teng et al. 2005), this is 
not much smaller than that found in radio-quiet QSOs (from —3 to —1). Moreover, absorp- 
tion may be a factor even at these energies. If the absorbing column exceeds ~10 24 cm -2 , the 
primary continuum emission is suppressed significantly by absorption and Compton down- 
scattering. Thus observations at > 10 keV are best to detect Compton-thick AGNs. The 
sensitivity of Suzaku at these high energies is well suited for the study of highly obscured 
sources like ULIRGs. 

In this paper, we present Suzaku XIS (0.5-10 keV) and HXD/PIN (15-40 keV) obser- 
vations of five well-known local ULIRGs. In §2, we discuss our sample. In §3, we report the 
observations and describe the methods we used to reduce the data. In §4, the results from 
our spectral analysis of the Suzaku data are discussed. In §5, we combine the Suzaku data 
with earlier published and unpublished XMM-Newton and Chandra data to fine tune our 
spectral models. The results of this study are summarized in §6. Throughout this paper, we 
adopt the cosmology of H =75 km s -1 Mpc -1 , f2M=0.3, and Q\=0.7. 



2. Sample 



The five ULIRGs in the present study are F05189-2524, F08572+3915, Mrk 273, 
PKS 1345+12, and Arp 220. They were selected because they have readily available Suzaku 
data, either from our own program (PI: Veilleux) or from the public archive. They are 
among the nearest, b rightest, and best-studied ULIRGs in the IRAS Bright Galaxy Survey 
( ISanders et al.ll2003l ). Table [1] lists the basic properties of these sources. Here we briefly 
review the relevant literature on each source. 
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F05 189-2524 is an unresolved late stage merger surrounded by tidal debris (Veilleux 
et al. 2002, 2006) with "warm" infrared colors. It is optically classified as a Seyfert 2, but 
near-infrared spectroscopy of this source reveals the presence of an obscured broad line re- 
gion ( BLR) at Paa (Veilleux et al. 1999a, b). Spectra from previous ASCA (IRisaliti et al. 



2000 



Severgnini et al.l 120011 ). XMM-Newton (jlmanishi fe Terashimal |2004j), and Chandra 



(IPtak et al.l 120031 ) observations flatten out above 2 keV and are best fit by an absorbed 
power law with A% ~ 0.5 — lxlO 23 cm -2 , T ~1.0— 1.9, and a thermal component with 
kT~0. 1-0.9 keV. The absorption-corrected 2-10 keV luminosity derived from the published 
XMM-Newton and Chandra data is ~10 43 ergs s . 

F08572+3915 is another "warm" ULIRG, con sisting; of a pair of interacting galaxies with 
nuclear separation of ~6 kpc (jVeilleux et al.ll2002l . and references therein). The northwestern 
nucleus is classified as a LINER. A 2-cm radio core coincides with the northwestern nucleus 



(INagar et al.ll2003l ). However, only an upper limit on its 2-10 keV flux exists in the literature 



7.6 x 10~ 12 ergs cm" 2 s" 1 - from the HEAO satellite flPolletta et al.lll996T ). 



Mrk 273 is in t he early phase of a merger where the nuclei are separated by 68 pc 
(jVeilleux et al.l 120021 ) and is optically classified as a Seyfert 2 (jVeilleux et al.l Il999aj ) . A 
[Si VI] 1.96 /xm featur e, a strong indicator of AGN activity, is detected in this "cool" ULIRG 
(jVeilleux et al.lll999bl ). Its radio flux falls above the radio-to-FIR correlati on of starbursts 



and a bright AGN-like radio core is detected on VLB A scale in this object (jLonsdale et al. 



19931 ) . The Chandra 0.5-2 keV X-ray spectrum of Mrk 273 is best explained with a MEKAL 
plasma with kT ~1.3 keV and its 2-10 keV spec trum is best fit by an absorbed power law with 



r ~1.0 and Ah ~10 23 cm" 2 (IPtak et all 120031 ) . The flat slope of the 2-10 keV c ontinuum 



may be a result of reflection, but the Chandra observation of IPtak et al.l (120031 ) shows a 
very weak Fe Ka emission line with an equivalent width of ~0. 09-0. 44 keV. Comparisons 
of observations by ASCA, BeppoSAX, and Chandra show that this source exhibits possible 
long term flux variability (jXia et al.l 120021 ) . 



PKS 1345+12 is yet ano ther "warm" ULIRG . It is in the early stage of a merger with two 
nuclei separated by 4.0 kpc (jVeilleux et al.l 120021 ) . Optically cl assified as a Seyfert 2 galaxy 
(jSanders et al.lll988l ; lKim et al.lll9981 ). inf rared observa t ions b y I Veilleux et al.l (jl997l ) suggest 
a buried BLR at Pact. Observations by lEvans et al.l (119991 ) have shown that the eastern 
nucleus has colors consistent with reddened starlight while the the western nucleus has 



extremely red colors ind icative of an optic al quasar. The western nuc 



with peak CO e mission (lEvans et al.l Il999h . a radio core (INagar et al. 



eus is also coincident 



2003), and 0.5-8 keV 



X-ray emission (jlmanishi fe Terashimal 12004 ) . According to llmanishi fe Terashimal (J2004J) , 
the Chandra continuum is consistent with that of an absorbed AGN (power law T=1.8 with 
iVn ~4.5xl0 22 cm -2 ). The Chandra data also showed a narrow Fe Ka emission line with an 
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equivalent width of ~0.13 keV. 



Finally, Arp 220 is by far the best-studied ULIRG due to its vicinity. Optically classified 
as a LINER, this "co ol" ULIRG is also an early merger with nuclei separated by ~0.4 kpc 
( jVeilleux et al.l 120021) . Multi- wavelength data suggest the presence of a black hole in the 
western nucleus (IDownes fc Eckartl 120071 . and references therein). Previous Chandra obser- 
vations detected both nuclei as well as e xtended soft X-ray emission from lobes and p lumes 
that extend beyond the optical galaxy ( Clements et al. 2002 ; McDowell et al.lbood ). The 
full Chandra band nuclear spectrum is best fit by a thermal MEKAL compone nt with kT 
0.8 keV and a flat power law with V ~1.1 absorbed by a colum n of ~10 21 cm -2 (IPtak et al. 



20031 ). Analysis of XMM-Newton data by llwasawa et al.l (120051 ) suggests the presence of an 
Fe K emission line with equivalent width of ~1.9 keV emanating from the western nucleus. 



3. Observations and Data Reduction 

The details of the Suzaku observations are listed in Table [2j Three of the sources 
(F05189-2524, Mrk 273, and PKS 1345+12) were of our own program (PI: Veilleux, Anabuki 
was co-PI for F05189-2524) and the other two, F08572 + 3915 (PI: Gallagher) and Arp 220 
(PI: the Suzaku Science Working Group, or SWG), were downloaded from the public archive. 
All of the observations were performed at the HXD aim point to increase the sensitivity of 
the HXD. The analysis of the data on F05189-2524, F08572+3915, Mrk 273, and Arp 220 
was performed with version 6.3.1 of HEASoft and CALDB version 20071016. PKS 1345+12 
was observed at the end of Cycle 2 and after the completion of the data analysis on the other 
objects in the sample. Thus, the analysis of the data on PKS 1345+12 was performed with 
the more up-to-date version 6.4 of HEASoft and CALDB version 20080401. 



3.1. XIS Data Reduction 

The XIS data reduction followed the guidelines provided in the Suzaku Data Reduction 
Guidd^- The data were screened following the version 2 data screening criteria. Due to 
nearby field sources, cleaned events for all objects except PKS 1345+12 were then extracted 
in circular regions with 1' radii (the minimum recommended region size) centered on the 
targets wherever possible. The PKS 1345+12 field is rather empty, so the extraction region 
has a radius of ~3.1'. Since the extraction regions are so large, Chandra ACIS data were 



http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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used as a check to ensure no other X-ray sources were included in the extraction regions. 
Background events were extracted in same-sized nearby source-free regions. The source 
spectra were binned to at least 50 counts bin -1 for PKS 1345+12 and at least 15 counts 
bin -1 for the others so that x 2 statistics can be used when modeling the spectra. 

The response (RMF) and auxiliary (ARF) files were produced using xisrmfgen and 
xissimarf gen. The ARF files were generated assuming 400,000 incident photons and the 
default grid spacing. 

In the modeling of the XIS data, the XIS1, XIS2, and XIS3 detectors are assumed to 
have the data processing version 2.0 cross-normalization factors of 1.065, 1.035, and 1.067 
with respect to the XIS0 detector, respective!}!^]. 



The HXD/PIN data were reduced following the guidelines provided by the Suzaku team. 
The HXD/PIN spectra for both sources were extracted after the selection of good time 
intervals (the "ANDed GTI" from both the data and the non- X-ray background provided 
by the Suzaku HXD team). Dead-time corrections (on the order of about 5%) were applied 
to the extracted source spectra. The extracted spectra were binned using grpphal^. 



The response files were provided by the HXD team. Due to changes in the bias volt- 
ages and the threshold over time, we used the first-epoch response file for observations of 
F05189-2524, F08572+3915, and Arp 220, second-epoch for Mrk 273, and fourth-epoch for 
PKS 1345+12. 

As its name implies, the non-X-ray backgrounds (NXB) from charged particles modeled 
by the HXD team do not include the cosmic X-ray background (CXB) which peaks within 
the energy range of the HXD /PIN. The CXB for each galaxy was modeled following a recipe^] 
provided by the HXD team. The simulated CXB is approximately 5% of the NXB. The NXB 
and CXB were added together using mathpha to provide total backgrounds for the HXD/PIN 
data. 



2 See: | ftp://legacy.gsfc.nasa.gov/suzaku/doc/xrt/suzakumem o-2007- ll.pdfl 

3 We binned the data using group 31 2 32 63 4 64 95 8 96 127 32 128 255 64. Each set of three numbers 
represent the channel range grouped and the number of channel bins in each group. This choice of grouping 
was used so that each bin contains approximately the same number of photons. 



3.2. HXD/PIN Data Reduction 




http: / /heasarc. gsfc.nasa.gov/ docs/suzaku/analysis/pinjcxb.html 
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At the time of this writing, the accuracy of the processing version 2.0 HXD/PIN back- 
ground model is 3.8%!E However, a ~10% offset in the version 2.0 background model was 
discovered for data taken between March and May of 2006 due to changes in the PIN ob- 
serving mode. For the affected data (F05189-2524 and F08572+3915), processing version 
1.2 HXD/PIN background were used as recommended by the Suzaku team. Dead time cor- 
rections were also performed on the backgrounds for the affected dataB The reproducibility 
of the version 1.2 background is between 5 and io%9. For the modeling of the data, the 
cross- normalization of the PIN with respect to XISO is 1.16 and 1.13 for the version 1.2 and 
2.0 backgrounds, respectively^. 



4. Suzaku Results 

The spectra were analyzed using the XSPEC package (version 11.3.2ag). The energy 
range of the detectors were limited to 0.5-10 keV for the XIS detectors and 15-40 keV for 
the PIN detector to avoid calibration problems. All the errors in the parameters are at the 
90% level for one parameter of interest (A% 2 = 2.7). 

In the modeling of the spectra, a simple power-law distribution absorbed only by the 
Galactic column is first assumed. If the model is not a satisfactory fit to the data, then 
more components are added on to the model until a satisfactory fit is achieved. The con- 
tinuum model components considered in this paper are: (1) an absorbed power-law model 
which represents emission from an AGN (model parameters fitted include a column den- 
sity, ]Vh, a power law index, T, and a normalization factor), (2) a thermal MEKAL model 
representing emission from stars or galactic winds (with a gas temperature, kT, and a nor- 
malization factor), (3) a partial-covering fraction model where some fraction, f cover , of the 
intrinsic nuclear radiation is absorbed and the rest goes through unimpeded. The other pa- 
rameters fitted are the equivalent column density, iVn, T, and a normalization factor. (4) a 
scattering model where an absorbed power law represents the direct transmitted component 
and an unabsorbed power law of the same photon index represents a scattered component 
(A^h, r, and a normalization factor for each of the direct and scattered components), (5) a 



5 This is the 1-cr statistical plus systematic error in the 15-40 keV band for a net integration time of 
10 ksec. See: | ftp://legacy.gsfc. nasa.gov/suzaku/doc/hxd/suz akumem o-2007-09.pdf. 



6 See: http: / /www. astro. isas.ac.jp/suzaku/ analysis /hxd/vl /pinnxb/ 

7 See: ftp://legacy.gsfc.nasa.gov/suzaku/doc/hxd/suzakumemo-2006-43.pdf 



8 See: [ftp: / /legacy.gsfc.nasa.gov/ suzaku/ doc/xrt / suzakume mo-2006-40.pdf for version 1.2 backgrounds 
and suzakumcmo-2007-ll.pdf for version 2.0 backgrounds. 
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neutral reflection modeQ with an unabsorbed power law representing the direct component 
and a reflection component with the same photon index as the direct component (T and a 
normalization factor for each of the reflected an d transmitted compone nts), and finally, (6) 
a pure ionized reflection model as proposed by iRoss fc Fabian without any direct 

transmitted component (T and a normalization factor). 



4.1. IRAS F05189-2524 

4-1.1. XIS Spectrum 

The XIS spectrum of F05 189-2524 is shown in Figure HJ First, we find that it is 
well-modeled by a simple power-law model plus a MEKAL component with abundances 
fixed at solar all modified by Galactic absorption (x„ = 1-57 for 124 degrees of freedom). An 
additional Gaussian component is needed (A% 2 =22.8 for a change in d.o.f. of 3) to reproduce 
the emission line at ~6.4 keV (i.e. the Fe Ka line; equivalent width, or EW, ~1.84 keV). 
While the spectrum shows excess emission at around 1.3 and 1.9 keV (Si XIII emission 
features), the addition of a second and a third Gaussian component to the model does not 
significantly improve the fit (Ax 2 =0.05 for a change in d.o.f. of 3). The excess emission 
at 1.9 keV may be a calibration feature associated with absorption edges from the mirror. 
The relatively flat spectral index and the large equivalent width of the iron line suggest that 
the spectrum may be reflection dominated. The addition of a reflection component to the 
continuum model results in a better fit (x 2 = 149 for 123 d.o.f.). The F-test probability for 
the addition of one parameter (the normalization of the reflected component) is relatively 
small, at 1.6 x 10 -2 , implying that it is reasonable to add the extra component. While the 
spectral index required for this fit, T = 2.68^013, is unusually steep for an AGN where the 
cano nical value for T is ~1 .8, it is not too different from the range seen in PG quasars (1.3- 



2.48; iPiconcelli et al.l 120051 ) . The best-fit reflection model parameters to the XIS spectrum 



are listed in Table [3] and the best-fit model is shown with the spectrum in Figure [H 



9 The reflection component is represented by the PEXRAV model in XSPEC. For the spectral fitting, the metal 
abundance, iron abundance, and inclination angle are fixed at the default value. The relative reflection 
parameter is fixed at -1, thus modeling only the reflection component. No energy cut off is assumed. 
Therefore, only T and the normalization factor are free parameters. 
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4.1.2. HXD/PIN Data 

F05189-2524 is undetected by the HXD/PIN. As mentioned in § 13.21 the documentation 
provided by the Suzaku team state that the reproducibility of the version 1.2 HXD/PIN 
backgrounds are on average ~5-10%. Here we estimate the systematic and statistical error 
of the HXD/PIN data for our specific observation. The standard deviation of the histogram 
of the 15-40 keV residuals from the 1-day (40 ksec) observations in the documentation imply 
that the systematic error of the NXB is ~2.4%. Considering the net HXD/PIN integration 
time and the 15-40 keV total background count rate for observations with integration times 
longer than 40 ksec, the statistical error of the F05189-2524 observation is ~0.9%. Thus, 
the total l-o" error (systematic and statistical errors added in quadrature) for the HXD/PIN 
observation of F05189-2524 is ~2.6%. Therefore, the 15-40 keV 3-a detection limit for this 
observation is ~0.020 counts per second (cps). 



4.2. IRAS F08572+3915 

4-2.1. XIS Non-Detection and XMM-Newton and Chandra Data 

As mentioned previously in § [21 there has never been a published X-ray detection of 
F08572+3915 in the X-ray energy range. It is also undetected by the XIS. A constraint 
on the limiting flux of F08572+3915 can be derived from the background. The average 
background count rate from the XIS front illuminated detectors (XIS0, 2, 3) is ~1.3xl0~ 3 
cps. Assuming an ideal, flat background (r ~1), the WebPIMMSo application was used to 
estimate the observed flux of the background spectrum. The application yields a flux of 
~4xl0 -14 ergs s _1 cm -2 , more than two orders of magnitude better than the upper limit 
placed on this object by HE AO. 

F08572+3915 was observed with both XMM-Newton and Chandra, but the results of 
the analysis were never published. F08572+3915 is undetected in the archived XMM-Newton 
data (ObsID: 0200630101, PI: Imanishi, t e fr e ctive=13 ksec). But the lower background and 
better spatial resolution of the Chandra ACIS-S archival data (ObsID: 6862, PI: Komossa, 
t e ffective=15 ksec) provide a conclusive detection of nine counts in the 0.5-8 keV band. The 
Chandra detection coincides with the optical northwestern nucleus (i.e. the LINER/Seyfert 2 
nucleus). Assuming a canonical AGN power law spectrum (r=1.8) modified by Galactic 
absorption for the Seyfert 2 nucleus, the detected count rate of F08572 + 3915 in the Chandra 



http : / /heasarc . gsfc . nasa.gov/Tools/w3pimms . html 
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band corresponds to an observed 0.5-10 keV flux of ~ 5 x 10 15 ergs s 1 cm 2 . This estimate 
is consistent with the upper limit derived from the XIS observations. 

The Chandra detection implies that the Seyfert nucleus may be a weak X-ray source 
or a heavily obscured AGN. If we assume no intrinsic absorption, the 0.5-10 keV flux from 
the Chandra observation of this object implies an intrinsic luminosity of ~5xl0 40 ergs s -1 . 
This luminosity falls withi n the range of X-ray luminosities of LINERs as measured by 
Terashima fc Wilson! (120031 ) and is consistent with LINERs being powered by low-luminosity 
AGNs. 



However, like many other LINER 2s (ITerashima fc Wilsonll2003f ). F08572+3915 is likely 
to be affected by absorption. Infrared observations of this object do show strong absorption 
features and weak PAH emission characteristic of deeply buried AGNs with a line-of-sight 
extinction of Ay > 78 mag. (Imanishi et al. 200 6; Armus et a l. 200 7; Veilleux et al. 2008). 
A CO-based estimate of the column density by lEvans et al.l (120021 ) is in the range of ~3- 
10xl0 24 cm -2 . The hardness rati o of the F08572+ 3915 detection is 0.56. Following the 
hardness ratio method presented in lTeng et al.l (120051 ). a single power law with an estimated 
photon index of ~-0.43 would fit the data. The inverted power law spectrum may be an 
indication of a large column density since the softer 0.5-2 keV photons are more readily 
affected by absorption. The estimated 0.5-10 keV flux from the hardness ratio method, 
which assumes Galactic absorption only, is 2.61 x 10~ 13 ergs s -1 cm -2 , corresponding to a 
luminosity of ~ 2 x 10 42 ergs s -1 . 



4.2.2. HXD/ PIN Data 

F08572+3915 is not detected by the HXD/PIN above the total background. Following 
the procedure outlined in § 14.1.21 the HXD /PIN 3-a detection limit was used to approximate 
the limiting count rate and observed flux in the 15-40 keV band. The 3-cr upper limit on 
the count rate for F08572+3915 in the HXD/PIN band is ~0.020 cps. 



4.3. Mrk 273 

4-3.1. XIS Spectrum 

Mrk 273 is detected by the XIS. We first tried to model the XIS spectrum with a 
MEKAL component describing the starburst and a heavily absorbed power-law distribution 
which does not fit the 2-10 keV spectrum properly. Rather, the XIS spectrum is well 
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modeled by the scattering model plus a MEKAL component with abundances fixed at solar 
and all components modified by Galactic absorption. Approximately 9% of the intrinsic AGN 
emission is scattered. An additional Gaussian component is needed to model the emission 
line at ~6.4 keV (EW~0.56 keV). While the spectrum shows excess emission at around 
1.9 keV (Si XIII emission), the addition of a second Gaussian component to the model does 
not significantly improve the fit (A% 2 =1 with Ad.o.f.=3). Again, the feature at 1.9 keV 
may be a calibration feature associated with absorption edges from the mirror. A reflection 
model was also tested on the XIS spectrum, but this results in a worse fit (xl — 1-52 for 187 
d.o.f.) than the scattering model (xl = 1-25 for 187 d.o.f.). The best-fit parameters to the 
XIS spectrum are listed in Table [3] and the XIS spectrum with the best-fit model is shown 
in Figure [H 



4.3.2. HXD/PIN Data & Contaminants 

Mrk 273 is marginally detected by the HXD/PIN. The net (observed minus background) 
HXD/PIN spectrum of Mrk 273 is ~5.4% of the total (NXB+CXB) background. Figure M 
is a comparison of the spectrum, total background, and net spectrum. The current best 
estimate for the error in the version 2.0 NXB background is 3.8% for observations with net 
integration time of 10 ksec. We estimate the error on our Mrk 273 observation following 
the steps outlined in §4.1.21 The total systematic plus statistical error for the Mrk 273 
observation is ~3.0%, implying the HXD/PIN detection of Mrk 273 is at ~1.8-cr (or with 
~93% confidence). 

Because the PIN is not an imaging detector and has a very large field of view (34' x 34', 
~3.5 times the area of the XIS field of view), there may be sources of contamination in the 
HXD/PIN signal. While hard X-ray and soft gamma ray catalogs (e.g. the INTEGRAL 
reference catalog; Ebisawa et al. 2003) do not list any sources of comparable flux within 
the PIN field of view and energy range, most of these flux values are based on modeling 
of previous observations at lower energies. There remains a possibility that the signal may 
come from a previously unobserved nearby source. Three likely candidates are detected 
within the field of view of the XIS. Mrk 273x is an unabsorbed Seyfert 2 galaxy at a redshift 
of 0.458 located approximately 1.2' northeast of Mrk 273. The XIS spectrum of Mrk 273x 
is well fit by a power-law component modified only by Galactic absorption (r ~1.69^q 
xl ~1.21 for 77 d.o.f.). The observed 0.5-10 keV flux of Mrk 273 is approximately 2.9 
times that of Mrk 273x. SBS 1342+560 is a background QSO at a redshift of 0.937 located 
approximately 6.3' south of Mrk 273. The XIS spectrum of this source is also well fit by 
a power-law component modified only by Galactic absorption (r ~2.03±0.05; xl ~l-28 for 
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242 d.o.f.). The observed 0.5-10 keV flux of Mrk 273 is approximately 1.5 times that of 
this source. Lastly, SDSS J1342512. 06+554759. 6 is a background QSO at a redshift of 1.17 
located approximately 7' southeast of Mrk 273. Its XIS spectrum is once again well fit by 
a power-law component modified only by Galactic absorption (T ~1. 62+0. 07; xl ~l-29 for 
139 d.o.f.). The observed 0.5-10 keV flux of Mrk 273 is approximately 1.9 times that of this 
source. 

Thus, none of the objects within the XIS field of view shows evidence of being highly 
obscured. One would therefore expect their spectral energy distributions at higher energies 
to follow the power law seen at low energies. If that is indeed the case, their contributions 
to the overall PIN signal of Mrk 273 are negligible; this is explained in the next section. 



4.3.3. Combined XIS-HXD/PIN Spectral Modeling 



We modeled the combined XIS and PIN spectrum of Mrk 273 with the scattering model. 
All of the best-fit parameters are consistent with the values obtained from the fit of the XIS 
spectrum alone with the exception of the column density which increased by a factor of ~2. 
The HXD/PIN contributions from the neighboring sources have been taken into account 
by adding the best-fit XIS model of the contaminants to the HXD/PIN component of the 
Mrk 273 model (see Figure [3D- The best-fit model to the full-band Mrk 273 spectrum in this 
case has parameter values listed in Table [3] and the spectrum is shown in Figure H] with the 
XIS-HXD/PIN spectrum. While the photon index value of ~1.4 may suggest flattening of the 
spectrum due to reflection, the addition of a reflection component resulted in a statistically 
worse fit (xl = 1-66 for 197 d.o.f.) relative to that obtained from the scattering model 
(xl = 1-35 for 199 d.o.f.). Approximately 5% of the intrinsic AGN flux i s scattered. The ga s 



temperature of ~0.7 keV is consistent with the range found in ULIRGs (jGrimes et al.ll2005l ). 



Franceschini et al. 


2003: 


Tene et al. 


2005. 


2008) 



(e.g. 

the XIS data with published results is postponed until § 15.21 



4.4. PKS 1345+12 

4-4- 1- XIS Spectrum 

The XIS spectrum of PKS 1345+12 is shown in Figure [H It is best modeled by a 
scattering continuum model (§ H]) modified by Galactic absorption. The transmitted AGN 
flux is ~2% of the intrinsic AGN flux. Unlike other ULIRGs, the 0.5-2 keV spectrum does 
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not require a MEKAL thermal component. The best-fit parameters are listed in Table [3j 
These best-fit values and t he model of the 0.5-10 keV con tinuu m are consistent with Chandra 
measurements reported by llmanishi fe Terashimal (120041 ) . The llmanishi fe Terashimal (120041 ) 
model for the Chandra data includes a narrow Fe line at 6.4 keV. This emission line also 
appears in the XIS spectrum, but is statistically insignificant (A% 2 =2.0 for a change in d.o.f. 
of 2). This is also consistent with the fact that the lower limit to the equivalent width of the 
Fe line in the llmanishi fe Terashimal (120041 ) model is keV. No flux or spectral variation were 
detected over the 0.5-10 keV energy range between the Chandra and Suzaku observations. 

The unabsorbed 0.5-10 keV flux for PKS 1345+12 is ~6.0xl0 43 ergs s" 1 and the ab- 
sorption corrected ratio log[L 2 _io kev/^m] is ~-2.3, consistent with the range found in PG 
quasars. This and the lack of a thermal component in the 0.5-2 keV spectrum suggest that 
the X-ray spectrum is dominated by an AGN. 



44.2. HXD/PIN Data 

The net HXD/PIN spectrum of PKS 1345+12 lies ~3.4% above the total background 
(Figure [2]). Accounting for the uncertainties in the PIN background of ~3.1% (see § 14.1.21 
on how this was calculated), the detection is at only l.ler. There are not enough net counts 
for a meaningful spectral fitting of the full-band (XIS+HXD/PIN) spectrum. However, 
assuming no contamination from nearby sources, the net observed 15-40 keV count rate for 
PKS 1345+12 is 9.9±2.8xl0~ 3 cps. Using the XIS best-fit model, given the lack of an obvious 
reflection component, the 15-40 keV count rate translates to a flux of ~4.4^3 x 10~ 12 ergs 

— 1 —2 

s cm . 



4.5. Arp 220 

4-5.1. XIS Spectrum 



Arp 220 is detected by the XIS. At first, the XIS spectrum was modeled by a simple 
power-law model plus a MEKAL component for the thermal contribution from the starburst. 
An additional Gaussian component is needed to reproduce the emission line at ~6.7 keV 
(EW~1.98 keV), which is consistent with emission from Fe XX-Fe XXVI. This model with 
only Galactic absorption is a satisfactory fit to the XIS data (xl ~1.03 for 174 d.o.f.). The 
XIS extraction window encompasses the two nuclei of Arp 220 and a significant fraction 
of the extended soft X-ray emission ([McDowell et al.ll2003l ). Nevertheless, the power-law 
plus MEKAL model to the XlS-only data are consistent with the results from the Chandra 
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observation of the nuclear binary source (jClements et al.l 120021 ). There is, therefore, no 
obvious sign of spectral or flux variability in Arp 220. 

The Suzaku data on Arp 220 can be interpreted in several ways: (1) an AGN is not 
present, (2) there is a low-luminosity, unobscured AGN, or (3) there is a heavily obscured 
AGN. 

In the first case where an AGN is not present, the detected Fe line at ~6.7 keV would 
originate from hot gas heated by the heavy starburst activity. However, the measured 
equivalent width for th e line of L98l^ ' ^ keV seems too large for an Fe emission line arising 
from a pure starburst. Ilwasawa et al.l (120051 ) fitted the 2.5-10 keV continuum of the XMM- 
Newton data with a collisionally ionized plasma model. The model requires a gas temperature 
of ~7 keV and a metallicity of ~2 times solar in order to reproduce the 6.7 keV Fe feature 
(EW~2 keV). This temperature is inconsistent with the detected Ca XIX line in the XMM- 
Newton data. The XIS data do suggest the presence of an emission feature at ~3.9 keV 
consistent with Ca XIX, but this detection is not statistically significant. 



The weak but unobscured AGN case is unlikely based on CO observations by lDownes fe Eckart 

1.3 x 10 25 cm" 2 for the 



( 120071 ). Their observations with IRAM imply a column density of 
western nucleus of Arp 220. Unless the covering fraction of the absorber is much less than 
unity, the unobscured AGN interpretation of the data is inconsistent with the CO observa- 
tions. 

In contrast, the third scenario of a heavily obscured AGN is consistent with the Fe 
line detection and th e CO observation of a large column density. The ionized reflection 
model as proposed by iRoss fe Fabian I (120051 . and references therein) can reconcile the large 
equivalent width of the highly ionized iron line as well as the large absorbing column. Since 
the absorbing column is so large, the observed spectrum is purely reflected. The XIS data 
is best fit by a MEKAL plus an ionized reflection model with the ionization parameter (£) 
fixed at 10 3 ergs cm s _1 (xl = 1-03 for 176 d.o.f.; see Tableland Figured]). Thus the third 
scenario for explaining the Suzaku data on Arp 220 is favored. 

Based solely on the XIS data, the unabsorbed 0.5-10 keV luminosity of Arp 220 is 
only ~ 1.8 x 10 41 ergs s _1 cm -2 , much lower than that found in quasars. The absorption- 
corrected ra|io_l_o£jX2-jn kev/^m] is ~-4.9, consistent with those found in some nearby 
ULIRGs (ITeng et al.ll2005l ). but lower than those of radio-quiet quasars. Of course, if Arp220 
is highly absorbed then the intrinsic log[L 2 _ 10 kov/^m] is higher and may be within the range 
for PG quasars. The thermal component of the spectrum has a 0.5-2 keV luminosity of 
~ 3 x 10 40 ergs s -1 . If this lu minosity is c o mplet ely due to thermal bremsstrahlung, then 



applying this to Equation (1) of ITeng et al.l (120081 ). the spatial extent of the emitting region 
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is ~ 2 — 11 kpc, depending on the choice of filling factor (10 3 < f gas < 10 x ) and assuming 
an electron density of 1 cm -3 . 

4.5.2. HXD/PIN Data 

As shown in Figure [2J the net HXD/PIN spectrum of Arp 220 is only ~2.2% above the 
(NXB+CXB) background , so it is within the uncertainties of the PIN background modeling 
(~3.0%, calculated following procedures in § I4.1.2p . Thus, the estimated 15-40 keV count 
rate for Arp 220 is ~0.025 cps. 

5. The Long-term Variability of IRAS F05189-2524 and Mrk 273 

Within the span of the Suzaku observations (~80 ksec each), no significant variability 
is detected in the XIS data for both F05189-2524 and Mrk 273. However, there is evidence 
for long-term variability in both sources when the Suzaku data are compared with previous 
ASCA, BeppoSax, Chandra, and XMM-Newton data. In this section, we discuss the results 
of these comparisons and their implications for the spectral models. 

5.1. IRAS F05189-2524 

F05189-2524 was observed by XMM-Newton in March 2001 and then by Chandra in 
October 2001 and January 2002. We extracted the spectra from the XMM-Newton and 
Chandra archives. There do not appear to be significant variations in the 0.5-2 keV flux of 
F05189-2524. The 2006 Suzaku measurement is within a standard deviation of the weighted 
average of the previous measurements. In contrast, the Suzaku 2-10 keV flux is a factor 
of ~30 lower than the measurements made by other observatories, ~4 standard deviations 
away from the weighted average of the previous measurements (see Figure [5] and [6]) . 

None of the XMM-Newton and Chandra data shows a significant iron feature near 
6.4 keV, in contrast to our more recent Suzaku data, so we did not consider reflection models 
for these archived data. The scattering model (§HJ) with a thermal MEKAL component gave 
adequate fits to the spectra. In the modeling of these data, we allowed the internal column 
density and the parameters for the power-law component describing the AGN within the 
source to vary freely. The models suggest that the archived XMM-Newton and the two sets 
of Chandra spectra are absorbed by a column density of 7.6 xlO 22 , 6.8 xlO 22 , and 6.1xl0 22 
cm -2 , respectively. We also deduce that the observed 2-10 keV flux of the AGN has de- 
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creased from 3x10 12 ergs s 1 cm 2 in the 2001 XMM-Newton data to 1x10 13 ergs s 1 
cm -2 in the 2006 Suzaku data i.e. a reduction of a factor of ~30. 

The decrease in flux and the change in spectral shape may be due to a number of 
reasons. Here, we test three scenarios by modeling the multiple epoch data simultaneously 
with a single model. The three scenarios considered are: (1) a change in the column density 
while the properties of the AGN remain the same, (2) a change in the covering fraction of 
the absorber while the properties of the AGN remain the same, and (3) the column density 
and covering fraction remain the same while the intrinsic AGN luminosity is changed. In 
this scenario, a reflection component is visible in the "low-state" AGN. 

The first scenario is tested by modeling all four data sets simultaneously by a single 
scattering model plus a MEKAL component modified only by Galactic absorption and a 
Gaussian describing the iron line. All of the model parameters are set to be the same for 
each of the data groups, except for the column density parameter which is allowed to vary 
independently. This results in a satisfactory fit to the data (% 2 = 1.26 for 503 d.o.f.). This 
model and the multi-epoch data are shown in Figure [7J and the best-fit parameters are listed 
in Table U 

For the second scenario, the multi-epoch data are simultaneously fit by a single partial 
covering fraction model. All of the model parameters are set to be the same between epochs, 
except for the covering fraction. The modeling of this scenario results in an unacceptable 
fit with x 2 ~ 8. The addition of another partial-covering absorber to this model results in 
a still unsatisfactory fit with xl — 2.2 for 498 degrees of freedom. Thus the partial-covering 
model is not a good description of the data. 

Lastly, the third scenario is examined. In this scenario, the active nucleus has switched 
off prior to the final set of observations, leaving behind a residual reflec tion component 



visibl e in the "low" sta te. This has been ob served previously in NGC 4051 (IGuainazzi et al. 



19981 ) and NGC 1365 (IRisaliti et al.ll2005l ). The model for this scenario is a combination 
of a MEKAL component, plus a reflection component on top of a scattering model. A 
Gaussian is also added to model the Fe line seen in the Suzaku data. For the Suzaku-epoch, 
the scattered component normalization is held fixed at zero. Similarly, for the non-Suzaku- 
epochs, the reflected component normalization is fixed at zero. Additionally, the absorption 
parameter in each of the non-Suzaku-epochs is allowed to vary independently while the 
Suzaku-epoch absorption parameter is held fixed at zero. Since the intrinsic AGN luminosity 
has changed, the normalization factors for the direct component are different for the Suzaku- 
and non-Suzaku-epochs. The modeling of the data with this scenario is rather successful 
with xl — 1-30 for 502 degrees of freedom. The best-fit parameters to this model are listed 
in Table H] and the model shown in Figure [7J 
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Both the increase in column density and the decrease in intrinsic luminosity of the AGN 
appear to be equally good descriptions of the data. However, in the first scenario, the A^ H 
value required to explain the data is > 2 x 10 24 cm -2 . At this limit, we can no longer model 
it as a pure absorption spectrum since scattering both in and out of the line of sight are 
becoming important. 

The thermal 0.5-2 keV luminosity is ~ 1 x 10 41 ergs s _1 . Following the methods pre- 
sented in § 14.5.11 the physical extent of the emitting region is ~4-17 kpc). The absorption- 
corrected 0.5-10 keV luminosity of the AGN detected in F05189-2524 at its " high" state 



(2000 -2002 observations) is ~3 x 10 ergs s , almost as luminous as some quasars (lElvis et al 



1994] ). The absorption-corrected ratio log[L 2 -io kev/^m] is ~~2.5, consistent with that found 



in radio-quiet PG quasars. The corresponding numbers for the "low" state as observed 
by Suzaku are ~8xl0 41 ergs s -1 and -4.1, respectively, well below that of radio-quiet PG 
quasars. 



5.2. Mrk 273 



Xia et all (12002 ) compiled a list of observed fluxes of Mrk 273 and Mrk 273x from ASCA, 
BeppoSax, and Chandra dating fro m 1994 to 200C j^j We add to this list the XMM-Newton 
observation from 2002 analyzed by iBalestra et all (120051 ) and our Suzaku observation from 
2006. Figure [8] shows the observed fluxes of Mrk 273 and Mrk 273x over the 12-year period. 
Since ASCA and BeppoSax did not have the spatial resolution to separate Mrk 273 from 
Mrk 273x, the flux points from 1994 and 1998 represent the sum of the two galaxies. As 
the figure shows, the total flux between Mrk 273 and Mrk 273x appears to have increased in 
2006 in the 0.5-2 keV band and appears to be in a "high" flux state in 2000 relative to all 
the other measurements. 

Although Mrk 273x was once th ought to be a B L Lacertae object because of its high 
X-ray to optical B-band flux ratio (IXia et all Il998l ). its 2-10 ke V flux has n ever s hown 
dramatic variability. The Suzaku data support the assessment of Kia et all (120021 ) that 
Mrk 273x is unlikely to be a BL Lac object. As discussed in §4.3.2, its spectrum is well fit 
by a simple powe r law m odified only by Galactic absorption, in agreement with the findings 



m 



Balestra et all (1200a ). 



Assuming a constant flux for Mrk 273x implies that the nominal observed 2-10 keV flux 



11 Although the authors tabulated the fluxes from 1996 to 2000, th e ASCA observ ation date was mis- 
identified and the data set was actually taken in December of 1994 (see llwasawa II1999T ) . 
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from Mrk 273 dropped by about a factor of two from 1994 to 1998, then increased by a factor 
of ~2.5 from 1998 to 2000, and then again dropped by a factor of more than two from 2000 
to 2002. It appears to have remained roughly the same from 2002 to 2006. This variability in 
flux is not very significant when considering the uncertainties of the measurements. However, 
the spectral variability is clearly seen when comparing the spectra of Mrk 273 from Chandra, 
XMM-Newton, and Suzaku (Figure [9]). The higher 0.5-2 keV flux in 2006 and the higher 
2-10 keV flux in 2000 as shown in Figure [8] correlate with the difference in the spectral 
shapes of the source in each of these epochs. 



5.2.1. Modeling the Multiple Epoch Data 

As with the multiple data sets of F05189-2524, we first modeled the archived Chandra 
and XMM-Newton spectra of Mrk 273 individually. These data sets were each best fit by the 
scattering model. The fits to the Chandra {XMM-Newton) data imply an intrinsic column 
density of 4.1 (6.8) xlO 23 cm -2 and an intrinsic 2-10 keV flux of 2.5 (2.0) xl0~ 12 ergs s~ 2 
cm -2 . A comparison of the Chandra [XMM-Newton) data with the fits to the 2006 Suzaku 
XIS data in Table [3] suggests that the intrinsic column density has increased by ~110% 
(~27%) while the intrinsic 2-10 keV luminosity has decreased by ~80% (0%) if we were 
to compare data only over the 0.5-10 keV range. However, the addition of the HXD/PIN 
data (Table [3]) implies that the column density has increased by a factor of ~4 (~2) and the 
intrinsic 2-10 keV luminosity has increased by a factor of ~3 (~4). 

Next, all three data sets are modeled with a single model simultaneously to determine 
whether the change in the spectral shape is due to the change in column density or the 
intrinsic AGN luminosity. Three scenarios are considered: the change in spectral shape is 
due to (1) a change in the absorbing column and the intrinsic AGN flux remains the same, (2) 
a change in the intrinsic AGN luminosity with the absorbing column remaining constant, and 
(3) a change in the covering fraction of the absorber with the column density and intrinsic 
luminosity remaining constant. 

To evaluate the first scenario, all three data sets are simultaneously modeled by a scat- 
tering plus MEKAL model modified only by Galactic absorption and a Gaussian describing 
the iron line. All of the model parameters are set to be the same for all of the data groups, 
but the column density for data sets from each epoch is allowed to vary (assuming the in- 
trinsic luminosity remains the same). As before, the contaminants in the HXD/PIN field 
of view are taken into account in the modeling. This interpretation of the data results in a 
poor fit (xl = 1-76 for 380 d.o.f.), much worse than the modeling to each data set alone and 
underestimates the flux of the HXD/PIN data (see Figure [TDT) . 
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The second scenario is explored by again modeling the multi-epoch data simultaneously 
with a single model. The model is the same as that used in the first scenario. However, 
in this case, the model parameters for each epoch are set to be the same except for the 
normalization of the power law. In this case, we are testing whether the change in the 
intensity of the power-law component can describe the observations. The best-fit model for 
this scenario takes into account the contributions from the nearby contaminants. While the 
model describes the 0.5-10 keV data for each epoch well (xt = 1-53 for 378 d.o.f.), it fits the 
HXD/PIN data very poorly, severely under-estimating the HXD/PIN flux (see Figure [T01 . 

To test the third scenario, the Chandra, XMM-Newton, and Suzaku spectra are modeled 
using two absorbers with the same column densities for data from each epoch but with 
covering fractions that are allowed to vary freely. The contributions to the HXD/PIN flux 
from the neighboring XIS sources are taken into account in the fitting of the Suzaku data. 
The best-fit model (xt = 1-39 for 376 d.o.f.) is shown in Figure [TU1 and the best-fit parameter 
values are listed in Table HI Based on Figure [10J the third scenario is the preferred model 
for the multiple epoch data. 

Using the parameters derived for this third scenario, the absorption-corrected luminosity 
from the AGN is 4.0 x 10 43 ergs s -1 in the 0.5-10 keV band and 2.3 x 10 43 ergs s -1 in the 
15-40 keV band. The ratio log[L 2 -io kov/^m] is ~-2.3, within the range found in radio-quiet 
PG quasars. The thermal 0.5-2 keV luminosity is ~ 1 x 10 40 ergs s _1 . Following the methods 
discussed in § 14.5-H the emission region is ~4-17 kpc in size. 

Given the lack of knowledge of the detailed geometry of the nuclear region, it is difficult 
to derive a quantitative measure of the distance from the central source to the absorber. 
However, the year-to-year time variability suggests that the absorber is close (~1 pc) to the 
nucleus. 



5.2.2. Comparison with OSSE Data 



Dermer et al.1 (119971 ) observed Mrk 273 with the OSSE instrument on-board the Comp- 
ton Gamma Ray Observatory. Though their observations resulted in a non-detection, they 
placed upper limits on the 50-100 keV b and flux and the column density. Possible expla- 
nations suggested by iDermer et al.1 (119971 ) for the non-detection include the possibility that 
this ULIRG is gamma-ray weak or the AGN source is highly variable in gamma-rays and 
they happened to have observed the object in the low-flux state. 

Another explanation for the non-detection is that the AGN in Mrk 273 is hidden behind 
a large and patchy column of gas > 10 24 cm -2 . This is not inconsistent with the upper limit 
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derived from previous CO measurements (< 2 x 10 24 cm -2 , assuming a covering fraction 
of unity; see Dermer et al. 1997 for more detail). These limits on the column density are 
consistent with our partial covering absorption model (JVh ~ 1.6 x 10 24 cm -2 ). From an 
extrapolation of the HXD/PIN observation and model, we derive an upper limit to the 50- 
100 keV photon flux of 6.2 x 10~ 5 photons cm -2 or 6.9 x 10~ 12 ergs s -1 cm -2 . This value is 
more stringent by a factor of ~2 than the OSSE value (1.5 x 10~ 4 photons cm -2 ). 

6. Summary 

The results of our analysis of Suzaku XIS (0.5-10 keV) and HXD/PIN (15-40 keV) obser- 
vations of five of the brightest and best-known local ULIRGs (F05 189-2524, F08572+3915, 
Mrk 273, PKS 1345+12, and Arp 220) have been presented and compared with earlier Chan- 
dra and XMM- Newton data. The results can be summarized as follows: 

1. The XIS observations of F05189-2524 reveal a significant change in the observed 2- 
10 keV spectrum relative to previous Chandra and XMM-Newton observations. The 
spectral variation in F05189-2524 suggests that the central source may have turned 
off, leaving behind a residual reflection component. However, an increase in column 
density cannot be completely ruled out. The absorption-corrected 2-10 keV luminosity 
to infrared luminosity ratios of F05189-2524 during its "high" state is consistent with 
values observed in PG quasars. 

2. The XIS spectrum of Mrk 273 contains a strong Fe Ka line and shows a change in the 
observed 2-10 keV spectrum relative to previous Chandra and XMM-Newton obser- 
vations. Mrk 273 is marginally (1.8-er) detected at high energies, with the HXD/PIN 
spectrum ~5.4% above the background. A change in the covering fraction of the ab- 
sorber best explains the spectral variations in Mrk 273, although a drop in the intrinsic 
AGN luminosity cannot be formally ruled out. A column density of ~10 24 cm" 2 is de- 
rived from the Suzaku data. The changes in spectral shape and covering fraction on a 
time scale of a few years suggest that the absorbing matter is ~1 pc from the central 
source. The Mrk 273 spectrum is best modeled by a ~94% covering fraction model. 
The absorption-corrected 2-10 keV luminosity to infrared luminosity ratios of Mrk 273 
is consistent with values observed in PG quasars. 

3. F08572+3915 is undetected in the Suzaku XIS and HXD/PIN observations. The low 
X-ray count rate derived from unpublished Chandra observations, combined with mid- 
infrared observations, suggests that this source is highly obscured. 
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4. PKS 1345+12 is detected by the XIS. No apparent flux or spectral variability is detected 
in its 0.5-10 keV spectrum relative to previous Chandra observations. The net 15- 
40 keV HXD/PIN spectrum is only l.lcr above the total background, not strong enough 
for a meaningful full-band (XIS+HXD/PIN) spectral fitting. Unlike other ULIRGs, 
the 0.5-2 keV spectrum of this source does not contain an obvious thermal MEKAL 
component. Combining this result and the fact that the absorption-corrected 2-10 keV 
to infrared luminosity ratio of PKS 1345+12 is in agreement with those of PG quasars, 
the data suggest that the X-ray luminosity of this object is dominated by an AGN. 

5. Arp 220 is detected by the XIS, but not by the HXD/PIN. Its 0.5-10 keV spectrum, 
including the iron complex, appears unchanged since previous Chandra observations. 
The X-ray continuum emission is in agreement with the possibility that a highly ob- 
scured AGN is present. This interpretation of the data is consistent with previous CO 
observations. The measurements of the iron emission line at 6.7 keV can be reconciled 
by the ionized reflection model. 

6. In all three cases where an optically thin thermal component is contributing to the soft 
X-ray emission detected by XIS (F05 189-2524, Mrk 273, Arp 220), the temperature, 
kT ~0.3-0.8 keV, is consistent with previous observations of ULIRGs. 
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Table 1. The Sample 



Object 


z 


log Fir 


^25 i-im/ f&0 fim 


Spectral 


Nh, Galactic 


Scale 


D L 






[crgs/s/cm 2 ] 




Type 


[10 20 cm" 2 ] 


[kpc/"] 


[Mpc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


F05189-2524 


0.042 


-8.87 


0.25 


S2/S1 


1.92 


0.71 


173.3 


F08572+3915 


0.058 


-9.18 


0.23 


LINER/S2 


2.60 


1.05 


242.0 


Mrk 273 


0.038 


-8.78 


0.10 


S2 


1.10 


0.70 


156.3 


PKS 1345+12 


0.122 


-9.65 


0.35 


S2/S1 


1.90 


2.05 


531.8 


Arp 220 


0.018 


-8.11 


0.08 


LINER 


4.27 


0.34 


72.9 



Note. — Col. (1): Object name. Col (2): redshift. Col. (3): logarithm of infrared (8-1000 
fim) flux. Col. (4): the 25-to-60 fim IRAS flux ratio, a measure of the du st temperature. Col. 
(5): galaxy spectral type. Col. (6): Galactic hydrogen column density llDickev fc Lockmanl 
1990). Col. (7): Physical size corresponding to 1". Col. (8): Luminosity distance. 



Table 2. Suzaku Observations 



Object 


Observation 


PI 


Date 


XIS Net Exposure HXD Net Exposure 




ID 


Name 


[UT] 


[ksec] 


[kscc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


F05189-2524 


701097010 


Veilleux; Anabuki 


2006 April 10 


78.2 


48.0 


F08572+3915 


701053010 


Gallagher 


2006 April 14 


77.2 


58.1 


Mrk 273 


701050010 


Veilleux 


2006 July 7 


79.9 


76.3 


PKS 1345+12 


702053010 


Veilleux 


2008 January 7 


53.0 


41.4 


Arp 220 


700006010 


Suzaku SWG 


2006 January 7 


98.6 


86.9 



Note. — Col. (1): Object name. Col (2)— (3): Suzaku proposal number and principal investigator of 
archived data. Col. (4): beginning observing date in UT. Col. (5)-(6): net exposure time in kiloseconds 
after screening and HXD deadtime corrections. 



Table 3. Best-fit Parameters to Suzaku Spectra^ 



Object 

(i) 



Best-fit 
Model 
(2) 



kT 
[keV] 

(3) 



(4) 



H. source 
> 23 cm" 2 ] 
(5) 



/ 

(6) 



E linc 

[keV] 
(7) 



[keV] 
(8) 



EW 
[keV] 

(9) 



(10) 



X /d.o.f. Observed Unabsorbcd Observed Unabsorbcd Observed Unabsorbcd 

F 0.5-2 keV F 0.5-2 keV F 2-10 keV F 2-10 keV F 12-40 keV F 12-40 keV 
(11) (12) (13) (14) (15) (16) (17) 



XIS Only 



F05189-2524 
Mrk 273 

PKS 1345+12 
Arp 220 



71+ - 13 
64 +0 04 

67+ 07 
U D -0.04 



2 ' 68 -0.13 
1 56+°' 15 

1 64+ 09 
1 '°^-0.11 

1SM -0.07 



1-71 
1.20 
, + 0.03 
' -0.07 

07+° 02 
uu '-0.02 



82.73 6A0±°H O.Oolg;^ 



8.643 
0.43 4 



0.09 6.39 
0.02 
— 6.70 



0.87 
+ 0.03 
0.03 



+ 0.13 
0.08 



-.+0.10 
J -0.00 



, + 0.27 
5 -0.05 



' -0.41 
g + 0.15 
°-0.11 



,+ 0.54 
-0.32 



2.20 
3.65 
3.13 
1.80 



+ 0.48 
0.74 
5+1.84 
0.90 
+ 0.52 
0.64 
+ 0.78 , 
0.40 ' 



X 10 

■ in-9 i 



185/123 
233/187 
164/156 
181/176 



0.66 
1.17 
0.37 
0.66 



+ 0.10 
0.10 
+ 0.10 
-0.13 
+ 0.09 
0.09 
+ 0.10 
0.11 



70 +0 ' 14 
'"-0.10 

a qk+ 410 

8 ' M8 -2.05 
5,93 -3.73 

-, or + 0.75 

L -' sb -0.i0 



1 03+ ' 39 
1UJ _0.38 

4 29+ 13 

9 53+ 317 
y - oo -4.06 

98+ 24 
u ' ao -0.42 



I 03+ 040 
1 ' ucl -0.40 

20 44+ 9 ' 48 
' !u - 44 -4 .84 

II 91 + 3 ' 12 

1 04+ 072 
" -0.35 



Combined XIS+HXD/PIN Spectra 


Mrk 273 


B 


0.64 4 : 9 !! 4 i.4oj:°;i° i-*.09±i;|| 0.05 6.3S±° g ° l o.03± g 8 3 o.es 4 :"; 1 ^ e.wj; 2 ; 3 ? x 10- 4 268/199 i.is±g i4.5oj:|;| 2 isstm 43 - l9 i2 e 9 3 6 


' 6 - 59 ±35.5l 


K 97 +28.79 
°" 4 '-21.45 



Note. — t F08572+3915 is undetected by the XIS and is thus not included in this table. * Since the spectrum is pure reflection, the normalization factor is of the reflection component. Col. (1): 
Object name. Col. (2): favored model — (A) reflection plus power law model, (B) scattering model, (C) ionized reflection model. Col (3): temperature of the thermal starburst component from the 
MEKAL model. Col. (4): slope to the power-law model. Col. (5): absorption within the source. Col. (6): the ratio of the normalization parameter between the reflected (or scattered) component to 
the direct component. Col. (7): rest energy of the emission line. Col. (8): line width of the emission line. Col. (9): equivalent width of the emission line. Col. (10): normalization factor of the ^^cct 
component in photons cm -2 s _1 at 1 keV. Col. (11): x 2 fitting statistics per degree of freedom. Col. (12)-(17): flux in units of 10 — 13 ergs s — 1 cm — 2 . " 



Table 4. Best-fit Parameters to Spectra of F05 189-2524 and Mrk 273 from Multiple Epochs 



Telescope 

(i) 


kT 

[keV] 
(2) 


r 

(3) 


/ 

(4) 


N H sourco 1 
[10 23 cm- 2 ] 
(5) 


/cover 1 N H aourco 2 
[10 23 cm" 2 ] 
(6) (7) 


/cover 2 B line CT 

[keV] [keV] 
(8) (9) (10) 


EW K 
[keV] 

(11) (12) 


Observed 
F 0.5-2 keV 

(13) 


Unabsorbcd 

F 0.5-2 keV 
(14) 


Observed Unabsorbcd Observed 
F 2-10 keV F 2-10 keV F 15-40 keV 

(15) (16) (17) 


Unabsorbcd 

F 15-40 keV 
(18) 














F05189-2524: 


scattering model, change in 


column density (x 2 /d.o.f. 


= 640/503) 









XMM (2001) 0.23+° ^| 1.91 + °-^ 0.02 0.58±o'o| 

Chandra (2001) 0.23 1.91 0.02 0.90+p'°| 

Chandra (2002) 0.23 1.91 0.02 0.71 + °^ 

Suzaku (2006) 0.23 1.91 0.02 



6.37 
6.37 
6.37 
6.37 



+ 0.27 7 3+0.31 64 +0.98 1 42 +0.17 xl0 -3 1 .02 + ° \1 30.00±i 



0.73 
0.73 
0.73 



0.64 
0.64 
0.64 



n 77 +0.08 
u - ' ' -0.08 
n H7 +0.09 
0-87_ . 09 

70+ 10 



; 53 29.38t 5 4 .35 42.03±«;f g 



30.00 
30.00 
30.00 



95 33 
M '°'-5.79 

27.69+ 7 '^ 

1 45+ - 28 



42.03 
42.03 
42.03 



F05189-2524: reflection dominated, change in AGN luminosity (x /d.o.f. = 659/502) 



XMM (2001) 0.24+°'° 2 
Chandra (2001) 0.24 
Chandra (2002) 0.24 
Suzaku (2006) 0.24 



-0.13 
1.79 

1.79 

1.79 



02 
02 
0.02 
6.82 



0.59 
0.93 
0.72 



+ 0.05 
0.05 
+ 0.09 
-0.09 
+ 0.07 
0.06 



2 + 0.17 „ 00 +0.46 , qc + 3.64 
^-0.22 uuu -0.00 i MO -1.00 



i.27t°; 2 t * io- 

1.27 x 10 -3 
1.27 x 10 -3 
2 - 18 l :36 X ltr 



1.02 
0.77 
0.87 
0.69 



+ 0.12 
0.10 
+ 0.10 
-0.11 
+ 0.10 
0.10 
+ 0.09 
0.10 



27 02"*" 1 " 55 

27.02 

27.02 

75+° 10 
u - '5_ .io 



29.92; 3 ; 3 

2S - 9 °t3.03 
28.16" 



1.06 



+ 3.14 
3.82 
+ 0.33 
0.14 



42.80i u: ° 5 3 

42.80 

42.80 

1 06+ ' 36 
1 ' UD -0.11 



Mrk 273: change in covering fraction (x 2 /d.o.f. — 521/376) 



-to 

-J 

I — 



XMM(2002) 
Suzaku(2006) 



-0.03 
0.79 



,.,+0.04 
,o -0.02 
1.73 



,r, qi+0.86 „ Rn +0.01 
15.91_ 73 U.»0_ Q2 



,.91 0.94±° \ 
,.91 0.94 + °'°; 



3 00+ - 23 
o.uu_ 20 



96+ - 01 
U ' MD -0.01 

0.81±° ° 2 2 
0.69 + „'° 2 



6 40+ - 01 
D '* u -0.13 

6.40 
6.40 



00+ ' 06 
0.00 



0.25+0.04 

0.25 
0.25 



, 24 +0.20 , -3 
0.07 ^'^-0.11 x lu 
2.24 x 10~ 3 
2.24 x 10~ 3 



„ 7 ,, + 2.92 , 7 „ q +4.25 

U ''3-0.15 4 ''2 J -2.26 

0. 85±o'2? 47.29 

1. U + °°l 47.29 



7 49+ 3 - 06 

4 -o9t?;l 7 

4.36 + ?' 3 ° 



s4.io±I;f 8 

84.10 
84.10 



Note. — Col. (1): Spacecraft used and year of the observation. Col (2): temperature of the thermal starburst component from the MEKAL model. Col. (3): slope to the power-law model. Col. (4): the ratio 
of the normalization parameters between the reflected (or scattered) component to the direct component. Col. (5)— (8): absorption within the source and corresponding covering fraction for the covering fraction 
absorber model. Col. (9): rest energy of the emission line. Col. (10): line width of the emission line. Col. (11): equivalent width of the emission line. Col. (12): normalization factor of the direct component in 
photons cm" 2 s — 1 . Col. (13)-(18): flux in units of 10" 13 ergs s" 1 cm -2 . 
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Fig. 1.— Spectra of F05189-2524 (top left), Mrk 273 (top right), PKS 1345+12 (bottom left) 
and Arp 220 (bottom right) with their respective best-fit models for all four XIS detectors. 
F08572 + 3915 is not detected by the XIS. The data are binned to at least 50 counts bin" 1 
for PKS 1345+12 and at least 15 counts bin -1 for the others. The horizontal axis is energy 
in the observer's frame. In F05189-2524, Mrk 273, and Arp 220, an emission line is detected 
near 6.4-6.7 keV, consistent with emission arising from neutral or ionized iron. The same 
emission line also appears to be present in PKS 1345+12, but the detection is not statistically 
significant (Ax 2 =2.0 for Ad.o.f=2). Other than PKS 1345+12, the 0.5-2 keV spectrum for 
each object has a thermal component with temperatures of ~0.7-0.8 keV, consistent with 
earlier results on ULIRGs. The best-fit parameters are listed in Table [31 
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Fig. 2. — HXD/PIN spectrum (open circles), together with the background spectrum 
(NXB+CXB; small crosses) and the net spectrum (open squares) of Mrk 273, PKS 1345+12, 
and Arp 220. The net spectrum for Mrk 273, PKS 1345+12, and Arp 220 is 5.4%, 3.4%, 
and 2.2% above the background, respectively. Mrk 273 is the only source that is detected 
above the background beyond the uncertainties of the background modeling (at the 1.8-0" 
level). The spectra of F05189-2524 and F08572+3915 do not lie above the background, so 
they are not shown here. 
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Fig. 3.— The net background-subtracted 15-40 keV HXD/PIN spectrum of Mrk 273 com- 
pared with the spectral models of possible contaminants in the PIN field of view. The solid 
line represents the best-fit model from the XIS-HXD/PIN spectral fitting and is a sum of 
the flux from Mrk 273 and the contaminants. The dashed lines represent the best-fit XIS 
models of the contaminants extrapolated to the HXD/PIN energy range. The contributions 
from these contaminants to the overall PIN signals are negligible. 



31 
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Fig. 4. — The combined XIS-HXD/PIN spectrum of Mrk 273 with the best-fit scattering 
model. The horizontal axis is energy in the observer's frame. The directly transmitted AGN 
flux is ~5% of the intrinsic flux. The iron line at ~6.4 keV is detected and a MEKAL 
model with gas temperature ~0.7 keV is needed to reproduce the soft X-rays. The cross- 
normalization of the HXD/PIN with respect to XISO is assumed to be 1.13. The best-fit 
parameters for the full band (XIS+HXD/PIN) modeling are listed in Table [3j 
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Fig. 5.— Observed 0.5-2 keV (top) and 2-10 keV (bottom) fluxes of F05189-2524 from 
1995 to 2006 as determined from ASCA, BeppoSAX, XMM-Newton, Chandra, and Suzaku 
data. While there has been little change in the 0.5-2 keV flux, the nominal 2-10 keV flux of 
F05189-2524 has decreased by a factor of ~30 since previous observations. The dashed line 
represents the average "high" state flux (measurements made prior to Suzaku) as weighted 
by the measurement errors while the dotted lines denote one standard deviation away from 
the mean. The 2-10 keV flux value as measured by Suzaku is ~4 standard deviations away 
from the weighted mea n. The ASCA (1995) and BeppoSAX (1999) values are drawn from 



Severgnini et al.l (120011 ); the XMM-Newton (2001) and Chandra (2002) values are drawn from 



our modeling of the archived spectra with the scattering model in this work (see §5.11) . The 
error bars for the 2006 value are within the data point. 



-33- 



0.0005 



■ft 



0.0001 r 



> 

(D 5xl0 -5 
w 



s 

o 

\ 

c 

o 



10 



-5 _ 



o 5xl0" 6 
X! 



io- 6 

5xl0" 7 



A O 



0.5 



.(St: ^ 



□ 



o □ 



o D "^o 

O OCb 
□ 




a March 2001 
□ October 2001 
O January 2002 
■ April 2006 



10 



channel energy (keV) 



Fig. 6.— Comparison of the unfolded spectra of F05189-2524 from the 2001 March XMM- 
Newton (open triangles), 2001 October Chandra (open squares), 2002 January Chandra (open 
circles), and 2006 April Suzaku (closed squares) observations. The unfolded spectra were 
created using the best-fit models to each individual spectrum. Below ~1.3 keV, there appears 
to have been negligible change in flux or spectral shape between the different observations. 
However, above ~1.3 keV, the spectral change is obvious. The iron line is prominent in 
the Suzaku data, but is not noticeable in the other observations. The best-fit models to the 
spectrum changed from a scattering-dominated scenario in 2001 and 2002 to a reflection- 
dominated scenario in 2006. This may be an indication that the central source has faded 
prior to 2006 or variations in the column density. 
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Fig. 7. — The two models to the Chandra- ACIS, XMM- Newton-EPIC pn, and Suzaku- 
XIS/PIN data used to explain the change in the 2-10 keV spectral shape of F05189-2524: 
(top) change in the absorbing column alters the spectral shape of the source, and (bottom) 
the AGN has switched off and has left behind a residual reflection component (see §5.11 for 
details). The response from each detector is folded in with the data. The horizontal axis 
is energy in the observer's frame. The cross-normalization factor is assumed to be unity 
for XMM-Newton and Chandra with respect to Suzaku-XISO. Both models have comparable 
reduced \ 2 values (1.27 for 503 d.o.f. versus 1.30 for 502 d.o.f.) and appear to fit the overall 
0.5-10 keV spectrum well. Neither model can be ruled out by the data. Table H] lists the 
results from this modeling. 
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Fig. 8. — Observed 0.5-2 keV (top) and 2-10 keV (bottom) fluxes of Mrk 273 (crosses) and 
Mrk 273x (open circles) from 1994 to 2006 as determined from ASCA, BeppoSAX, XMM- 
Newton, Chandra, and Suzaku data. Each asterisk represents the sum of the fluxes from 
Mrk 273 and Mrk 273x. In the top panel, the dashed line represents the weighted average 
of the data points for Mrk 273 and Mrk 273x from 1994, 2000, and 2002 and the dotted 
lines represent a single standard deviation away from that mean. In the bottom panel, the 
dashed line represents the weighted average of the data points for Mrk 273 and Mrk 273x 
from 1994, 1998, 2002, and 2006 and the dotted lines denote one standard deviation away 
from the mean. Small flux increases in Mrk 273 may be observed in 2006 in the 0.5-2 keV 
band and in 20 00 in the 2-10 k eV band. The ASCA (1 994) and BeppoSAX (1998) values 
are drawn from llwasawa I ( 1999 ) and Risaliti et al. ( 2000 ). respectively. The Chandra (2000) 
and XMM-Newton (2002) values are derived from our modeling of the archived spectra with 
the scattering model (see §5.2p . The error bars to the total flux are the sums of the errors 
in the measured fluxes of Mrk 273 and Mrk 273x added in quadrature. A small horizontal 
offset was applied to the Mrk 273 data points to better display the measurement errors. 
Since the errors on the BeppoSAX measurements are not published, they are assumed to be 
±20% based on the results on F05189 - 2524 bv ISevergnini et all J200lh . 



-36- 



0.0005 



> 
X 

VI 

B 

o 

\ 

C 

o 
o 

£1 

Oh 



0.0001 
5xl0- 5 



: a A A 

o 



io- 



5x10 



-6 _ 



a - B -£.%r 




o 



o 


2000 


A 


2002 


■ 


2006 



10 -6 LL- L 



A o® 



o 




■s 



a t)*- A . .j5<# - . 



A AD 



O 



A ■ 

A A 



A A 
A" 



-i i u 



Zfi" ■ ?? 
■ 

A 

_l I I l_ 



0.5 



1 5 
channel energy (keV) 



10 



Fig. 9. — Comparison of the unfolded spectra of Mrk 273 from the 2000 Chandra (open 
circles), 2002 XMM-Newton (open triangles), and 2006 Suzaku (closed squares) observations. 
The unfolded spectra were created using the best-fit models to each individual spectrum. 
The shapes of the spectra appear consistent between 1 and 2 keV. The flux of the Suzaku 
spectrum is higher than that of the others below 1 keV, while the Chandra spectrum appears 
to be higher than the others between 4 and 6 keV. These differences are correlated with the 
flux variability as shown in Figure This comparison also shows that the iron line at 
6.4 keV is detected in all three observations. While the 2-10 keV flux variability is not 
very significant, the spectral shape of the source has changed between 2000 and 2006. This 
change may be due to the variations in the column density. 
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Fig. 10. — The three models to the Chandra- ACIS, XMAf-Newton-EPIC pn, and Suzaku- 
XIS/PIN data used to explain the changes in the 2-10 keV spectral shape of Mrk 273: 
scenario 1 (top) tests only a change in the absorbing column, scenario 2 (middle) tests only 
a change in the intrinsic luminosity of the AGN, and scenario 3 (bottom) tests the change in 
the covering fraction of the absorber (see §5.21 for details). The response from each detector 
is folded in with the data. The horizontal axis is energy in the observer's frame. The cross- 
normalization factor is assumed to be unity for XMM-Newton and Chandra with respect 
to Suzaku-XISO. While all three models have similar reduced \ 2 values (1.76, 1.53, and 
1.39 respectively), the third model provides the best-fit to the overall 0.5-40 keV spectrum. 
Thus, we favor the third scenario (change in the covering fraction of the absorbers) as the 
best explanation for the change in the 2-10 keV spectral shape of Mrk 273. Table H] lists the 
results from this modeling. 



